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Genetic susceptibility is involved in the pathogenesis of vitiligo. Association studies with a 
whole genome-based approach instead of a single or a few candidate genes may be useful 
for discovering new susceptible genes. Although the etiology of non-segmental and 
segmental types is different, the association between gene polymorphisms and vitiligo has 
been reported, without defining types or in non-segmental type. Whole genome-based 
single nucleotide polymorphisms (SNPs) were examined in patients with non-segmental 
and segmental types of vitiligo using the Affymetrix GeneChip 500K mapping array, and 
10 functional classes of significant SNPs were selected. Genotyping and data analysis of 
selected 10 SNPs was performed using real-time PCR. Genotype and allele frequencies 
were significantly different between both types of vitiligo and three of the target SNPs, 
DA//4/75(rs2277046), SJRN3 (rs2273171), and KIAA 1005 (rs321 3758). A stronger 
association was suggested between the mutation in /C/A4?005(rs3213758) and the 
segmental type compared to the non-segmental type of vitiligo. DNAH5 (rs2277046), 
STRN3 (rs227317l), and KIAA1005 (rs3213758) may be new vitiligo-related SNPs in Korean 
patients, either non-segmental or segmental type. 
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INTRODUCTION 

Vitiligo is an acquired depigmentary disorder of the skin. Al- 
though the etiology is unclear, the melanocytes involved in viti- 
ligo are genetically susceptible and vulnerable to extracellular 
insults. Identification of genes involved in susceptibility to vitili- 
go has been attempted by association studies of candidate genes, 
which related to pathogenesis of non-segmental vitiligo includ- 
ing autoimmunity (1-8). However, genetics of vitiligo is charac- 
terized by incomplete penetrance, multiple susceptibility loci 
and genetic heterogeneity. Therefore, association studies with a 
whole genome-based approach instead of single or a few can- 
didate genes may be useful to discover new susceptible genes. 
With completion of the human genome sequencing, there is a 
focus on developing new sequencing methodologies consider- 
ing read length, thorough-put, and cost (9). A technical improve- 
ment in whole genome amplification with small amount of ge- 
nomic DNA has been adopted for single nucleotide polymor- 
phism (SNP) analysis (10). Recently, genome-wide association 
analysis has been attempted in vitiligo, with identification of 
genes including X-box binding protein 1 (XBP1) (11), SMOC2 
(12), and TYR (13). A study could be conducted with a unique 
isolated population with elevated prevalence of generalized vit- 
iligo and other autoimmune disease, which may provide a con- 



dition for detection of susceptible alleles with smaller sample 
sizes (12). However, studies have frequently been done with in- 
dependent patients and controls (11, 13). 

There is a consensus that vitiligo should be classified into 
segmental and non-segmental vitiligo forms based on clinical 
manifestations (14). The etiology of each type is considered to 
be distinct, although the association between segmental and 
non-segmental vitiligo has been reported. Nevertheless, the as- 
sociation between gene polymorphisms and vitiligo has been 
reported without defining types. In cases in which the type has 
been clarified, the association has usually been identified in 
non-segmental form. 

A genome-wide association study was conducted in Korean 
patients with vitiligo, although the whole genome SNP analysis 
was performed in a very small number of patients. Functional 
classes of significant SNPs were validated by genotyping pa- 
tients and independent healthy controls. Association between 
SNPs and vitiligo types was examined. 

MATERIALS AND METHODS 
Subjects 

Twenty patients with non-segmental vitiligo were examined for 
whole genome-based SNPs. In total, 163 Korean patients with 
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vitiligo were included to validate the 10 selected target SNPs in 
this study. A total of 1 13 cases were non-segmental, and 50 were 
segmental. An additional 97 patients with non-segmental (71 
patients) and segmental vitiligo (26 patients) were added, re- 
sulting in 184 patients with the non-segmental and 76 with the 
segmental type to analyze the association between the three 
significant gene mutations and either type of vitiligo. In the 184 
non-segmental and 76 segmental types, 86 and 27 patients were 
male, and 98 and 49 were female, respectively. Their ages rang- 
ed from 5 to 81 yr (mean, 40.1 yr) and 4 to 70 yr (mean, 23.6 yr), 
respectively. The control group included 192 healthy individu- 
als (82 males and 110 females), with ages between 30 and 80 yr 
(mean, 50.0 yr). 

Genome-wide genotyping with significant SNP analysis 

Genomic DNA was isolated from peripheral blood leukocytes 
using the GenEX Genomic kit (Geneall, Seoul, Korea). Whole 
genome-based SNPs were examined using 2 ug of the genomic 
DNA from 20 patients with non-segmental vitiligo. The Affyme- 
trix GeneChip 500K mapping array, which contained 500,568 
probes, was used for the genome-wide association analysis. The 
genotypes were determined using the BRLMM algorithm. Be- 
fore applying the algorithm, whole genome-based SNPs from 
44 individuals offered by the Affymetrix were put together for 
the genotype calling to compensate for the small number of pa- 
tients. The rates of mismatch between the genotyping from 20 
patients and that from 20 patients plus the 44 individuals ranged 
from 0.14 to 0.39%. The genotype calls from 20 patients with vit- 
iligo were compared with SNPs from 192 Korean healthy con- 
trols. All SNPs with a call rate less than 0.5, missing genotype 
frequency more than 0.5, and minor allele frequency less than 
0.01 were removed with any monomorphic SNP. Any SNP with 
a P value less than 0.0001 in Hardy- Weinberg equilibrium test 
was also removed to analyze valid SNPs. Risk alleles were ana- 
lyzed by the chi-square test using dominant, recessive, co-dom- 
inant, and additive genetic models. Significant SNPs were cho- 
sen by the level of significance which was either a = 0.001 or a = 
0.01 with multiple test correction, a Bonferroni method. 

Target SNP selection 

The number of significant SNPs was 1 14. They were divided into 
eight functional classes, such as non-synonymous, synonymous, 
mRNA UTR, intron, locus region, in-gene, intergenic, and un- 
known. The number of non-synonymous SNPs was eight, which 
were selected as the target SNPs. The SNPs within the linkage 
disequilibrium (LD) block, which included non-synonymous, 
synonymous or mRNA UTR SNPs, were also examined and 
seven significant SNPs were detected. Two (rsl7035120 and 
rs2588400) of the seven SNPs were already included as target 
SNPs. One SNP (rs3763979) with a function class of synonymous 
and the other (rsl444216) with a function class of in-gene were 



Table 1. Target SNPs selected by genome-wide association analysis 



No. 


Gene symbol 


Chromosome 


rs number 


Functional class 


fevel of 
significance 


1 


PLB1 


2 


2199619 


Non-synonymous 


a = 0.001* 


2 


DNAH^ 

Ui vrtl l\J 


5 


2277046 






3 


KIAA1005 


16 


3213758 






4 


STRN3 


14 


2273171 






5 


GLT8D2 


12 


17035120 






6 


L0C285556 


4 


13114983 






7 


L0C388559 


19 


2870128 






8 


L0C642624 


12 


2588400 






9 

10 


SLC16A7 
MUC19 


12 
12 


3763979 
1444216 


Synonymous 


01 = 0.01* 



The level of significance; The level of significance with multiple test correction, a 
Bonferroni method. 

added to the list of target SNPs, resulting in a total of 10 target 
SNPs (Table 1). 

Genotyping for clinical validation of the target SNPs 

The real time-polymerase chain reaction (PCR) using genomic 
DNA from patients with vitiligo and healthy controls was per- 
formed to examine the association between non-segmental or 
segmental vitiligo, the 10 target SNPs, and the three significant 
SNPs. 10 ng of each genomic DNA was mixed with 20 uL of re- 
action mixture containing 1.0 pL of reaction mix, 2.0 pL of Fast- 
Start DNA master, and 1.6 pL of MgCl 2 . Real-time PCR was done 
automatically using specific primers and probe included in the 
LightSNiP kit (TIB MOLBIOL, Berlin, Germany) with a LightCy- 
cler 480 Instrument (Roche, Mannheim, Germany). The Light- 
SNiP kit was designed as an individual SNP base by TIB MOL- 
BIOL. 

Statistical analysis 

All statistical tests were two-sided, and the level of significance 
was 0.05. Substitution of missing values was not performed. 
Data from both vitiligo and healthy control groups were com- 
pared using the chi-square test. One allele association analysis 
and four genotype association analyses were performed using 
dominant, recessive, co-dominant, and additive genetic mod- 
els. Odd ratios were calculated by chi-square test from a 2 x 2 
contingency table. All output was produced using SPSS version 
9.1 (SPSS Inc., Chicago, IL, USA). 

Ethics statement 

The institutional review board of Dongguk University Ilsan Hos- 
pital approved this study (IRB No. 2005-2), and the study was 
conducted according to the Declaration of Helsinki Principles. 
After obtaining informed written consent, blood was taken for 
isolation of genomic DNA. 

RESULTS 

The genotype frequencies of the 10 target SNPs were compared 
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Table 2. Genotype and allele frequencies of each gene polymorphism in the vitiligo patients and the healthy controls 



SNP 


benotype/Allele 


Vitiligo 
(n = 260) 


Normal 
(n = 192) 


Model 


Odd ratio (95% U) 


p* 


UNAn rszz/ /U4b 


A A 

AA 


en /oo cn\ 


67 (34 


90) 


U 


i QO M 0-\ O 77\ 

1 .00 (1 .Z\ -z./ /) 


U.UUo 




AT 
Mb 


1 1 Q MR 77\ 


89 (46 


35) 


p 
n 


9 nn ft oq q 1 o\ 

Z.UU { I .Zo-j. I L.) 


n 009 
u.uuz 




nn 


9,0 CW R^ 

OZ (O 1 . JM-y 


36 (18 


75) 




9 RQ /1 R^ A Xf\ 


n nnn 
u.uuu 




A 


9^7 MR Rftt 


223 (58 


07) 


AH 
ML] 


1 Rfl /1 9^ 9 C\8\ 

I ,DU \ I C-.\J0) 


n nnn 
u.uuu 




u 


9£^ fR4 49^ 


161 (41 


93) 


A 


1 RR /1 97 9 1 fi^ 

I ,DJ ( I .C-t I 


n nnn 
u.uuu 


CTDMO rnOOTQi 

olHNJ rszz/oi l\ 


TT 
1 1 


G1 /OO 

b 1 (zo.4b) 


59 (30 


73) 


n 
U 


1 a c /n nc; o on\ 
1 .40 (U.yo-Z.ZU) 


U.Uo4 




1 L- 


1 /rh nn\ 


103(53 


.65) 


D 

n 


1 .yo ( 1 .z 1 -0. 1 4j 


n nns 
u.uub 




cc 


69 (26.54) 


30 (15 


63) 


c 


2.22 (1.27-3.89) 


0.005 




T 


252 (48.46) 


221 (57 


.55) 


Ad 


1.47 (1.11-1.93) 


0.006 




C 


268 (51.54) 


163 (42 


.45) 


A 


1.44 (1.11-1.88) 


0.007 


rtZ*M7005rs3213758 


GG 


84(32.31) 


101 (52 


.60) 


D 


2.33 (1.58-3.42) 


1.69 x 10- 5 




AG 


77 (29.62) 


81 (42 


19) 


R 


11.19(5.65-22.18) 


4.48 x 10 12 




AA 


99 (38.08) 


10(5.21) 


c 


11.90 (5.84-24.25) 


9.06x10 12 




G 


245(47.12) 


283 (73 


70) 


Ad 


2.54 (1 .94-3.32) 


8.88 x 10 12 




A 


275 (52.88) 


101 (26 


30) 


A 


3.15(2.37-4.18) 


3.31 x 10 15 



*Pvalue derived from a chi-square test. D, dominant model; R, recessive model; C, co-dominant model; Ad, additive model; A, Allele analysis. 



Table 3. Genotype and allele frequencies of each gene polymorphism in the vitiligo patients with non-segmental (NSV) or segmental (SV) type and the healthy controls (N) 



SNP 


Geno type/ 


NSV 


SV 


Model - 


Odd ratio (95% CI) 


p. 




Allele 


(n = 1 84) 


(n = 76) 


N vs NSV 


NvsSV 


N vs NSV 


NvsSV 


D/V/Wrs2277046 


AA 


44 (23.91) 


15(19.74) 


D 


1.71 (1.09-2.68) 


2.18(1.15-4.13) 


0.020 


0.017 




AG 


88 (47.83) 


31 (40.79) 


R 


1.71 (1.05-2.77) 


2.83 (1 .57-5.08) 


0.030 


0.001 




GG 


52 (28.26) 


30 (39.47) 


C 


2.20 (1.24-3.89) 


3.72 (1 .78-7.80) 


0.007 


0.001 




A 


176(47.83) 


61 (40.13) 


Ad 


1.48 (1.12-1.97) 


1.96 (1 .35-2.86) 


0.007 


0.001 




G 


192(52.17) 


91 (59.87) 


A 


1.51 (1.13-2.02) 


2.07 (1.41-3.03) 


0.005 


0.000 


S77?/V3rs2273171 


TT 


44 (23.91) 


17(22.37) 


D 


1.41 (0.89-2.23) 


1.54 (0.83-2.86) 


0.139 


0.173 




TC 


92 (50.00) 


38 (50.00) 


R 


1.91 (1.14-3.17) 


2.06 (1 .09-3.89) 


0.013 


0.026 




cc 


48 (26.09) 


21 (27.63) 


C 


2.15(1.18-3.91) 


2.43 (1.12-5.28) 


0.013 


0.025 




T 


180(48.91) 


72 (47.37) 


Ad 


1.44 (1.07-1.94) 


1.55 (1.05-2.31) 


0.016 


0.029 




C 


188(51.09) 


80 (52.63) 


A 


1.42 (1.06-1.89) 


1.51 (1.03-2.20) 


0.018 


0.033 


KIAA 1005 rs321 3758 


GG 


66 (35.87) 


18(23.68) 


D 


1.98 (1.31-3.00) 


3.58 (1.96-6.52) 


0.001 


3.2 x 10" 5 




AG 


53 (28.80) 


24 (31 .58) 


R 


9.94(4.91-20.11) 


14.73 (6.75-32.17) 


1.6 x 10- 10 


1.5 x 10' 1 




AA 


65 (35.33) 


34 (44.74) 


C 


9.95 (4.77-20.73) 


19.08 (8.03-45.32) 


8.8 x 10"™ 


2.4 x 10 ' 1 




G 


185(50.27) 


60 (39.47) 


Ad 


2.34 (1.75-3.13) 


3.97 (2.60-6.06) 


8.4x10- g 


1.9 x 10 10 




A 


183(49.73) 


92 (60.53) 


A 


2.77 (2.04-3.76) 


4.30 (2.89-6.39) 


6.2 x 10~ 11 


5.9 x 10 " 



*Pvalue derived from a chi-square test. D, dominant model; R, recessive model; C, co-dominant model; Ad, additive model; A, Allele analysis. 



between vitiligo patients and normal healthy controls. The gen- 
otype and allele frequencies showed a significant difference in 
three of the target SNPs, such as DNAH5 (rs2277046), STRN3 
(rs2273171), and KIAA1005 (rs3213758) in all four different ge- 
netic models (Table 2), regardless of the vitiligo type difference 
(Table 3), although the STRN3 genotype was not significantly 
different in the dominant model. In DNAH5 (dynein axonemal 
heavy chain5), the genotype and allele frequencies of 12658 
A > G (4220 Thr > Ala) in patients with vitiligo were significantly 
different from those in healthy controls. The AA genotype fre- 
quency was lower, whereas GG was higher in vitiligo patients, 
suggesting that G allele at the DNAH A12658G polymorphism 
increased susceptibility to vitiligo. The genotype and allele fre- 
quencies of STRN3 (striatin calmodulin binding protein 3) 1376 
T > C (387 Asn > Ser) were also significantly different between 
the patients and healthy controls. The TT genotype was less fre- 
quent, whereas CC was more common in patients than con- 



trols, indicating that the C allele at the STRN3 T1376C may in- 
crease susceptibility to vitiligo. In KIAA1005, the genotype and 
allele frequencies of 3854 G > A (1264 Asp > Asn) in vitiligo pa- 
tients were significantly different compared to healthy controls. 
The GG frequency was lower, whereas AA frequency was higher 
in vitiligo, showing an A allele at the KIAA1005 G3854A may in- 
crease susceptibility to vitiligo. 

The genotype and, particularly, allele frequencies of DNAH5 
12658 A > G and STRN3 1376 T > C were not significantly dif- 
ferent between the non-segmental and segmental forms in all 
four models (Table 3), whereas those of KIAA1005 3854 G > A 
was different (Table 3). Therefore, a direct comparison between 
non-segmental and segmental vitiligo was performed. In fact, 
the frequency of allele A in the KIAA1 005 gene was significantly 
(P = 0.025) higher in the segmental compared to that in non- 
segmental vitiligo, although the genotype frequency suggested 
only a tendency for statistical significance in recessive (P = 0.058), 
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Table 4. Genotype and allele frequencies of each gene polymorphism between non-segmental (NSV) and segmental (SV) 



SNP 


benotype/Allele 


(N 


NSV 
= 184) 


SV 
(N = 76) 


Model 


Odd ratio (95% CI) 


P* 


UNAn rszz/ /U4b 


A A 

AA 


A A 

44 


(23 


91) 


1 c; /m ~7A\ 
lo (ly./4) 


U 


0.60 


(0 


34 


-1.06) 


U.U/o 




AO 
Ab 


RR 
oo 


(47 


83) 


o I (4U. / y) 


D 

n 


0.78 


(0 


40 


-1.51) 


U.4DO 




nn 

uu 


Oil 


(28 


26) 




r 


0.61 


(0 


33 


-1.12) 


0 119 

U. I I tL 




/\ 


1 76 


(47 


83) 


R1 (AO 1 31 


AH 


0.75 


(0 


51- 


■1 .08) 


0 124 




G 


1 92 


(52 


17) 


01 *q Q7\ 




0.73 


(0 


50- 


-1.07) 


0110 




TT 
1 1 


44 


(23 


91) 


01 /OO 07\ 


n 
U 


0.92 


(0 


.51- 


-1 .69) 


u./yo 




TP 

1 u 


0,9 

yz 


(50 


00) 


oo pU.UU) 


D 

n 


0.92 


(0 


48 


-1 .73) 


u. / oy 




cc 


48 


(26 


09) 


17(27.63) 


c 


0.88 


(0 


41 


-1.89) 


0.748 




T 


180 


(48 


91) 


80 (52.63) 


Ad 


0.94 


(0 


64 


-1.37) 


0.749 




C 


188 


(51 


09) 


72 (47.37) 


A 


0.94 


(0 


64 


-1.37) 


0.749 


/OW005rs3213758 


GG 


66 


(35 


87) 


18(23.68) 


D 


0.67 


(0 


39 


-1.16) 


0.156 




AG 


53 


(28 


80) 


24 (31 .58) 


R 


0.55 


(0 


30 


■1 .02) 


0.058 




AA 


65 


(35 


33) 


34 (44.74) 


C 


0.52 


(0 


27 


■1.01) 


0.055 




G 


185 


(50 


27) 


60 (39.47) 


Ad 


0.73 


(0 


53 


■1.01) 


0.060 




A 


183 


(49 


73) 


92 (60.53) 


A 


0.65 


(0 


44 


-0.95) 


0.025 



"Pvalue derived from a chi-square test. D, dominant model; R, recessive model; C, co-dominant model; Ad, additive model; A, Allele analysis. 



co-dominant {P = 0.055), and additive {P = 0.060) models (Ta- 
ble^. 

DISCUSSION 

This genome-wide association study in Korean patients with vit- 
iligo identified that three new SNPs, such as DNAH5 (rs2277046), 
STRN3 (rs2273171), and KIAA1005 (rs3213758), were associat- 
ed with increased susceptibility to both non-segmental and 
segmental types of vitiligo (Table 2, 3), although the number of 
patients was too small to identify the actual number of signifi- 
cant SNPs. 

Although an association with other variants in DNAH5, STRN3, 
and KIAA1005 with particular diseases has been reported (15- 
21), a clinical association of the SNPs, such as DNAH5 rs2277046, 
STRN3 rs2273171, and KIAA1005 rs3213758, has not been re- 
vealed. The DNAH5 gene is located in the coding axon of chro- 
mosome 5pl5.2. The A12658G (rs2277046) variant had an ami- 
no acid threonine residue (Thr 4220 ) changed to an alanine resi- 
due (Ala 4220 ). Mutation in the DNAH5 is found in patients with 
primary ciliary dyskinesia (PCD) (15), a rare disease transmit- 
ted as an autosomal recessive trait and characterized by recur- 
rent airway infections due to abnormal ciliary structure and 
function. Primary defects in the structure and function of sen- 
sory and motile cilia result in multiple ciliopathies (16). Multi- 
ple variants of the KIAA1005 gene, which is located in coding 
exon of chromosome 16ql2.2 and called as a retinitis pigmen- 
tosa GTPase regulator-interacting protein 1-like {RPGRIP1L) 
gene, have also been associated with certain clinical manifesta- 
tions, particularly ciliopathies as in DNAH5 with neurological, 
renal and ocular manifestations (17). The D1264N (rs3213758) 
mutation had the amino acid aspartic acid residue (Asp 1264 ) 
changed to asparagine (Asn 1264 ). Since cilia may be a sensitive 
and easily damaged organelle (18), performing diverse biologi- 



cal functions as hair-like structures extending from the cell mem- 
brane, certain defects of the cilia in either melanocytes or kera- 
tinocytes may relate to vitiligo development, although further 
studies are required. STRN3 gene locates in the coding exon of 
chromosome 14ql2. The STRN3 rs2273171 variant induces the 
amino acid asparagine residue (Asn 387 ) to change to a serine 
residue (Ser 387 ). Striatin is an intracellular protein characterized 
by four protein-protein interaction domains, a caveolin-bind- 
ing domain, a coiled-coil structure, a Ca 2+ -calmodulin (CaM)- 
binding domain, and a WD repeat domain (19). Results includ- 
ing the defective calcium uptake in both keratinocytes and me- 
lanocytes (20) and the activated nitric oxide synthases by in- 
creased levels of intracellular Ca 2+ via calmodulin (21), suggest 
an association of striatin with vitiligo development. 

In general, segmental vitiligo has been classified separately 
from other forms of vitiligo including non-segmental vitiligo. 
An association between gene polymorphisms and vitiligo has 
been reported for non-segmental vitiligo (1-4). In this study, the 
association was separately analyzed both in non-segmental 
and in segmental forms. The association did not differ in either 
type of vitiligo (Table 3), although a stronger association be- 
tween KIAA1005 (rs3213758) and segmental vitiligo was sug- 
gested (Tables 3 and 4). Neural dysfunction in the pathogenesis 
of segmental vitiligo and KIAA genes for organization and func- 
tion in the brain (22) may support this association. However, 
increasing clinical and etiological evidence suggests an associa- 
tion between segmental vitiligo and non-segmental vitiligo. 

Collectively, our genome-wide association study in Korean 
patients with vitiligo resulted in three new vitiligo-related SNPs, 
suchasD/VAH5(rs2277046), S77«V3(rs2273171), andKIAA1005 
(rs3213758), although a stronger association between KIAA1 005 
(rs3213758) and segmental compared to the non-segmental type 
of vitiligo was suggested. Future functional studies of these SNPs 
may contribute to the understanding of vitiligo pathogenesis. 
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